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Abstract 

We evaluate the on shell form factors of the electron for arbitrary momentum transfer and finite 
electron mass, at two loops in QED, by integrating the corresponding dispersion relations, which involve 
the imaginary parts known since a long time. The infrared divergences are parameterized in terms of a 
fictitious small photon mass. The result is expressed in terms of Harmonic Polylogarithms of maximum 
weight 4. The expansions for small and large momentum transfer are also given. 
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1 Introduction 



The imaginary parts of the on shell form factors of the electron at two loops in QED (the relevant graphs 
are shown for completeness in Fig.l) were calculated analytically long ago PP in terms of Nielsen's poly- 
logarithms |2] of maximum weight w = 3. The real parts however are still missing, as the integration of 
the dispersion relations for the form factors could not be carried out in closed form within the family of 
Nielsen's polylogarithms only. 

We show in this paper that the previous calculation can be completed within the family of the Harmonic 
Polylogarithms (HPL's) introduced in PJ. The result involves HPL's of maximum weight w = 4. 

The plan of the paper is as follows. In Section 2 we recall, from £Q, the proper dispersion relations 
satisfied by the two form factors. In Section 3 we list the results of the analytic integration for arbitrary 
momentum transfer t — —Q 2 and finite electron mass m. For positive Q 2 all the terms are real; as the 
expression is exact, it can be continued to t timelike and above threshold, where the form factors develop 
their imaginary parts. In Section 4 we give the expansions of the form factors for Q 2 — + oo, and in Section 
5 the expansions for Q 2 — > 0. In the Appendix A we recall the definition and the main properties of the 
HPL's, while in Appendix B we derive the analytic formulas which were used in the calculation. 




(a) (b) (c) 




(d) (e) 
Figure 1: The Vertex Graphs at two loops in QED (multiplicities understood). 



2 Dispersion Relations 



(2) 

The imaginary parts of the two on-shell vertex form factors of the electron at two loop in QED, ImF{ (£), 

(2) m 

Imi 7 ^ (t) were evaluated in for arbitrary value of the momentum transfer t and finite electron mass, 
within the Pauli-Villars regularization scheme for the ultraviolet divergences (needed for the renormalization 
of the inserted one loop subgraphs) and by using a small fictitious mass A for the regularization of the 
infrared divergences. In the imaginary parts t is above the threshold 4m 2 , where m is the electron mass; 
by introducing the dimensionless variable x through 

2 (l+z) 2 V~t- Vt-im 2 
t = m , x = — — ; , (2.1) 
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all the imaginary parts are can be written in terms of x, logx, log(l — x), log(l + x) and Nielsen polyloga- 
rithms [2] of maximum weight 3 and of argument x, —x, x 2 . 

As discussed in Eq.(1.30) of PP, the corresponding real parts are then given by properly subtracted 
dispersion relations, which for spacelike momentum transfer t = —Q 2 read 

Q 4 f°° dt' t' -Am 2 1 ,2), , . 

' -ImFf'tt'), (2.2) 



Q 2 + Am 2 J 4m 2 t> + Q 2 t' 2 ir 



Q 2 r dt' t'-Am 2 1 T (2) 



-lmF^(t'). (2.3) 



Q2 + 4m 2 y 4m2 t , + Q2 t , n 

Indeed, the above imaginary parts are singular at t = Am 2 , due to the infrared singularities of the 2-photon 
ladder graph showing up in the A — ► limit - hence the subtraction at Q 2 = —Am 2 and the factor (t' — Am 2 ) 

for making the dispersive integrals convergent at threshold. Note that the charge slope F± (0) and the 

(2) 

magnetic anomaly F 2 (0) are needed for writing the above dispersion relations (which therefore cannot be 

used, in that form, for their calculation). Let us further observe that Ea. l|2.2|l gives F± (0) — 0, the usual 
renormalization condition of the charge form factor. 
Similarly to Ea. l|2.1|l . in the spacelike case we put 



,2 



[1-yf ^Q 2 + Am 2 



V ^JQ 2 + Am 2 



(2.4) 



with the change of integration variable t' = m 2 (l+x') 2 jx' , the dispersive integrals occurring in Ea.s l|2.2l2.3(l 
are seen to be of the form 



dt' 

4m 2 ^ ^ 



r f(t')= fdx' (±j--±r t^tTt) ^')- 

Jo x' + y x' + l/yj 



(2.5) 



When the actual explicit expressions of the imaginary parts in terms of Nielsen's polylogarithms are used 
within the dispersive integrals Eg. Q2. 212. 31) . it is found that the resulting integrals cannot be evaluated 
in terms of that same class of functions. The Nielsen's polylogarithms, however, are special cases of the 
Harmonic Polylogarithms HPL's; when the more general formalism of the HPL's is used, one finds that 
not only the imaginary parts, but also the real parts given by the dispersive integrals can all be expressed 
in terms of HPL's of argument y and 1 (the HPL's of argument 1 are expressed in turn as combinations of 
a few mathematical constants like C(3), the Riemann £- function of index 3, 7r 2 and In 2). 

We give in the following paragraphs the explicit results of the integration. The relevant formulae are 
shortly discussed in the Appendix. 



3 Results 

We present in this section the results of the analytic integration of Eq. (|2.2() and Eq. I|2.3|l as a function of 
Q 2 = —t. As already said, A is the "small photon mass" used for the parameterization of the infrared 
divergences. All the polylogarithms depend on Q 2 through the variable y defined in Ea. (|2.4|l . Let us recall 
that y(Q 2 = oo) = 0, y(Q 2 = 0) = 1; when the momentum transfer t is positive (timelike) and varies from 
to the physical threshold t = Am 2 , Q 2 = — (t + ie) varies from —it to Q 2 = —Am 2 — ie; correspondingly, 
y = e 1 ^ , with <fi = at t = and (j> — it at t — Am 2 (and the form factors are still real despite the 
complex value of y). When t is above the threshold, it is convenient to write y — —x + ie, where x is the 
variable defined in Ea. (|2.1|) . so that as t varies from 4m 2 to +oo, x varies correspondingly from 1 to 0. For 
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y = —x + ie the polylogarithms with rightmost index equal to develop an imaginary part which generates 
the imaginary parts of the form factors (see 00] and the Appendix A for more details) . 
The explicit results are 
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f 2 (2) (-Q 2 



(to) { 



A \ f m 2 r 

7 



+ 



m 



#(0,0; y) 
23 



Q 2 + 4m 2 



#(0,0; y) 



VQ 2 (Q 2 + 4m 2 ) 



#(0;y) 



Q 2 v/Q 2 (Q 2 +4to 2 ) 



jS0 7r 4 + ^(0,-l; ? y)-^ 7 r 2 ff(0,0;y) 

(0, -1, 0, 0; y) + H(0, 0, -1, 0; y) - |tf (0, 0, 0, 0; y) 



- §C(3) + j^tf (-i ; y) - ^ 2 # (0; y) - \h$, o ; y) + ^(-1, o, o ; y) 

-if-if (0, 0, 0; y) - H(0, 1, 0; y) - £T (1, 0, 0; y) 



(Q 2 + 4m 2 ) 2 V<9 2 (<2 2 +4m 2 ) 



19 



7r 4 -7^ 2 tf(0;y) + 84C(3)tf(0;y) 



11 



(Q 2 + 4m 2 ) 2 



+157r 2 ff(0,-l;y)- y7 r 2 tf(0,0;y) 

-127r 2 i?(l, 0; y) - 42if (0, 0, 0; y) 

-126F (0, -1, 0, 0; y) + 192tf(0, 0,-1, 0; y) 

-3i?(0, 0, 0, 0; y) - 48ff(0, 0, 1, 0; y) 

+180ff(0, 1, 0, 0; y) - 72ff(l, 0, 0, 0; y) 

^tt 2 - 42 C (3) - ^H(-l; y) \^H{Q; y) ^H(0, 0; y) 

33 

+63if (-1, 0, 0; y) - 96H(0, -1, 0; y) - ytf (0, 0, 0; y) 



+24ff(0,l,0;y)- 90tf(l, 0,0; y) 



+ 



m 



(Q 2 + 4m 2 ) v /Q 2 (Q 2 + 4m 2 ) 



-4. 2 + i^ 7 r 4 -|i/(0 ; y) + ^ 2 i/(0;y) 

-28C(3) J ff(0; y) - 96tf (-1, 0; y) - ^7r 2 if (0, -1; y) 
+60tf (0, 0; y) + y 7r 2 ff(0, 0; y) + 24tf (1, 0; y) 
+47r 2 i7 (1, 0; y) + 7tf (0, 0, 0; y) + y H(0, -1, 0, 0; y) 



5 



+ 



m 



Q 2 + Am 2 



-21if(0, 0,-1,0; y) - -ff(0, 0,0,0; y) 
-20F(0,l,0,0;y) + 24F(l,0,0,0;j/) 



17 1 7 17 47 

y - 2tt 2 In 2 + -7T 2 + -C(3) + - n 2 H (-1; y ) - ^^^(0; y) 

95 1 19 

+ (0, 0; y) - - A H{-1, 0, 0; y) + 3H(0, -1, 0; y) - — ff(0, 0, 0; y) 



+ 



^Q 2 (Q 2 +4m 2 ) 



-H(0,l,0;y)-H(l,0,Q;y) 



j^ 2 + 2 -^H(0; y) + %n 2 H(0; y) + *±H(-1, 0; y) 



57 I I 

-H(0, 0; y) - 8H(1, 0; y) + —H(0, 0, 0; y) 



(3.7) 



4 Behaviour for Q' 
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We present the behaviour of the two Form Factors for Q 2 — > oo (corresponding to y — > 0), writing for short 
L = ln(Q 2 /m 2 ). The analytic continuation to large time-like t can be carried out with the replacements 
Q 2 = -{t + ie), L = \n{t/m 2 ) - in - so that for instance L 2 = ln 2 (i/m 2 ) - 7r 2 - 2in \n(t / m 2 ) . 
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The above formula matches with the sum of Eq.s(2.23) and (2. 34-2. 35a) of [5], which corresponds to the 
real part of F^ 2 \t) for large timelike t expanded in m 2 /t up to the zeroth order. 
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5 Behaviour for Q 2 — > 

We present the behaviour of the two Form Factors around Q 2 = 0, corresponding to y = 1. 




(5.10) 



(5.11) 



From the above expansions one can easily recover some familiar results, such as the renormalization con- 
dition F^ 1 (0) = and the two loop values for the slope (0) (the coefficient of the first term in Q 2 in 
Ea. (|5.1U|l with an overall minus sign) and of the electron anomaly F^iO). 
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A The Harmonic Poly logarithms, HPL's 



1-x ' 

1 



We recall for convenience of the reader the definition of the HPL's [3]. The HPL's form a family of functions 
depending on an argument, say x, and on a set of indices, say dj, i = 1, .., w or a in more compact notation, 
where each of the at can take one of the three values 1, 0, —1 and whose number w is called the weight of 
the HPL. At weight w = 1 there are 3 HPL's, defined as 

H(l;x) = -ln(l - x) , 
H(0;x) = lux , 

= ln(l+x), (A-l) 

whose derivatives can be written as 

^H(a;x)=f(a;x) , (A-2) 
where the index a can take one of the three values (1, 0, —1) and the rational factors f(a; x) are given by 

f(i;x) = 1 
f(P;x) = x 

/(-i;*) = j^. (A-3) 

At weight w > 1, if all the u> indices are equal to let us indicate them by W and define correspondingly 

H{Q w ;x) = —\n w x- (A-4) 

in all the other cases {i.e. when the indices are not all equal to zero), let us indicate any set of w indices by 
(a, 6), where a can take one of the three values (1, 0, —1) and b stands for the set of the other w — 1 indices, 
and define correspondingly 

H{a, b;x)= [ dx'f(a; x')H(b; x') . (A-5) 
Jo 

Note that in full generality (i.e. also when all the indices are equal to 0) one has 

^-H (a, b- x) = /(a; x)H(b; x) , (A-6) 
ax 

which can also be written as the equivalent indefinite integration formula 

dx'f(a; x') H(b; x') = A+ H{a, b; x) , (A-7) 

where A is an integration constant. 

Further, the product of two HPL's of the a same argument x and weights p, q can be expressed as a 
combination of HPL's of argument x and weight r = p + q, according to the product identity 

H{p-x)H{q-x) = Y, H ( f ;x), (A-8) 

r—pWq 

where p, q stand for the p and q components of the indices of the two HPL's, while pttl g represents all 
possible mergers of p and q into the vector f with r components, in which the relative orders of the elements 
of p and q are preserved. The simplest cases of the above identities are 

H(a;x)H(b;x) = H(a,b;x)+ H(b,a;x) , 
H(a;x) H(b,c;x) = H(a,b,c;x)+ H(b,a,c;x)+ H(b,c,a;x) ; (A-9) 
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more complicated cases are immediately established recursively (all the above formulae can indeed easily 
be checked by differentiating, repeatedly when needed, with respect to x). All the 3 W HPL's of weight w 
are linearly independent; Ea. ljA-9|) can however be used for replacing an HPL of weight w with products of 
HPL's of lower weight (but such that the sum of their weights is equal to w) and other HPL's of weight w. 

In the analytic continuation from the argument y in the range (0, 1) to y = —x + ie, where x is again in 
the range (0, 1), one finds that any HPL with rightmost index equal to does develop an imaginary part. 
It is therefore convenient to exploit Ea. l|A-8|) for expressing those HPL's in terms of HPL's having either 
no 0's on the right or only 0's as indices. One has for instance 

H(-l, 0; y) = H(-l; y)H(0; y) - H(0, -1; y) , (A-10) 

and similar formulae for more general cases. From the very definition one has further H(0; —x + ie) = 
H(0;x)+in, and H(-l\ -x + ie) = H(-l; -x) = -H(l;x), H(0, -1; -x + ie) = -H(0, l;x) , (those last 
two functions have no imaginary part in that range of values of x), so that finally one obtains 

H(-l, 0; -x + ie) = -H(0; x)H(l; x) + H(0, 1; x) + inH(l; x) . (A-ll) 

For a more complete discussion see 00]. 

Of particular interest are the values of the HPL's of argument equal to 1. They can be expressed as 
combination with rational factors of a limited number of mathematical constants, such as for instance 
Riemann £- functions, it 2 and In 2. 



B The definite integrals occurring in the calculation 

The imaginary parts of the form factors at two loops contain Nielsen's polylogarithms of maximum weight 
3 and arguments x' , — x' , x' 2 , which can be expressed in terms of HPL's of argument x' . Indeed, one 
has Li 2 (x) = H(0,l;x), Li 2 (-x) = -H(0,-l;x), Li 3 (x) = H (0,0,1; x), Li 3 (-x) = -#(0,0, -l;x), 
Si 2 (x) = H(0, 1, 1; x), S 12 (-x) = H(0, -1, -1; x) and 

S"i 2 (x 2 ) = 2H(0, 1, 1; x) - 2H{0, 1,-1; x) - 2H(0, -1, 1; x) + 2^(0,-1,-1; x). 

After the insertion of the explicit analytic expressions of the imaginary parts in Ea.s (|2.2l2.3|> and full partial 
fractioning in the integration variable x' ' , the integrands are found to consist of powers of the rational factors 
x' , 1/(1 + x 1 ), 1/(1 — x'), which depend on x' only, and \ j(x' + y), l/(x' + 1/y), which depend as well 
on y, times HPL's of argument x' and maximum weight w = 3. 

The integrals with the factors depending on x' only can be evaluated by parts or by means of the formulas 
defining the HPL's; they give end-point values of HPL's of argument 1 and maximum weight w = 4 (some 
care may be needed in properly grouping terms whose end-point values are otherwise separately divergent). 

For the analytic evaluation of the integrals involving l/(x' + y) and l/(x' + 1/y), along the lines of 
Section 7 of JHJ, we introduce two families of related functions depending, as the HPL's, on a set of w 
indices and the argument y, defined for w = 1 as 

F(-l;y) := f ^—=H{-l; V )-H{Q;y) , (B-1) 

Jo x +y 

r 1 dr' 

G(-l;y) := / -— -=H(-l;y), (B-2) 
Jo + - 

and for w > 1 as 

F(-l,b;y) 
G(-l,b;y) 



1 dx' 



H(b;x') , 



o x' + y 
1 dx' 

^-rH(b-x'). 



o x' + i 

V 



(B-3) 
(B-4) 
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Note that the first index of the above functions is frozen to the value —1 (a convention suggested by the 
y = 1 limiting value). 

It is easy to see, proceeding by induction on w, that the functions F(— l,a;y) and G{— l,a;y) with a 
total of w indices are just homogeneous combination of HPL's functions of weight w of argument y and of 
their values at 1. 

For w = 1, the result is obvious from the definitions Ea.s (|B-llB-2|) . Next, assume that the identities 
are established for F(— 1, b; y) up to a certain weight w — 1 > 1 and consider the function of weight w 

P 1 dx' 

F(-l,a,b;y) = -?—H(a,b;x'). (B-5) 



Jo x' + y 

Its values at y = 0, 1, according to the definition of the HPL's Ea. l|A-5p are 

F(-l,a,b;l) = H{-l,a,b;l) , 

F(-l,a,b;Q) = H(0,a,b;l) . (B-6) 
By differentiating Ea. (|B-5jl with respect to y and integrating by parts on x' one obtains 



d f 1 dx' 

—F(-l,a,b;y) = f(-l;y)H(a,b;l)- f(a;x') H(b;x') . (B-7) 

°V Jo x i V 

(We have assumed H (a, b; 0) = 0, which is the most frequent case; for H(a, b; 0) ^ see the remarks below.) 
In the case a — 0, after partial fractioning in x' , one has 

|-F(-l,0,6;y) = f(-l;y)H(0,b;l)-f(0;y) [ dx' ( f (0; x') - H(b;x') 
ay Jo V x' + y J 

= f(-l;y)H(0,b;l)-f(0;y)(H(Q,b;l)-F(-l,b;y)], (B-8) 



where Ea. (|A-5|l and Ea. ljB-3|) have been used to carry out the x'-integration. 
For a — 1,-1 one obtains similarly 

£-F(-l,l,%y) = -f(-l;y)F(-l,b;y) , (B-9) 
dy 

^F(-l,-l,b;y) = f(-liy)H(-l,b;l)+f(V,y)(H(-l,b l l)-F(-l,b l y)\ (B-10) 

In the r.h.s. of the previous equations appears the function F(—l,b;y) of weight w; when proceeding by 
induction in w one can substitute the already obtained identities expressing it in terms of iJ's of weight 
w — 1 and argument y. A final quadrature in y, carried out according to Ea. ljA-5|l . gives F(—l,a,b;y) in 
terms of HPL's of weight w up to an additive constant, which can be fixed by one of the Ea.s (|B-6l . 

It is not difficult, following the above lines, to work out all the formulas needed for the calculation. As 
an example, we give one of the identities of weight w = 4, 

F(-l, 1, 0, 0; x) = +-H(-1; x)H(0; x)H(Q; x)H(0; x) + 2H(-1; x)H(0,-l; 1)H(0; x) 
6 

-H(-l;x)H(0, 0,-1; x) + H(-l; 1)H(0, 0,-1; 1) + ff(-l; 1)H{0, 0, 1; 1) 
-iff (0, -1; x)H(0; x)H(0; x) + hi (0, -1; x)H(0, -1; x) 

-2H(0, -1; l)H(0, -1; x) + ^(0,-1; 1)^(0, -1; 1) + H(0, 0, -1; x)H(Q; x) 
-iJ(0,0,0,-l;x) + ff(0,0,0,-l;l)-iJ(0,0, 1,-1; 1) . (B-ll) 

Essentially the same procedure applies to the other family of functions G{— 1, a; x) as well. 

For particular values of the arguments, the end-point values occurring in the derivation may give diver- 
gent contributions when taken separately; in those cases, one can parameterize the end-point singularities 
by integrating on x' from S to 1 — i] and then take the 5, 7/ — > limit in the final result (which is of course 
finite) . 
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